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Abstract 

Masses of the lowest multiplets of c, 6-mesons in S- and P-wave states 
with quantum numbers J^*^ : 1 ,0+^,1"' ,1+^,2++ are ob- 
tained with the help of dispersion N/D method of heavy quark effec- 
tive theory. The results are in a good agreement with experimental 
data. Radiative decay widths for some states mentioned above were 
calculated. 

PACS numbers: 11.55.Fv, 12.39.Ki, 14.40.Lb, 14.40.Nd 
1. Introduction 



The mass spectrum of the charmed and beauty mesons is the subject 
of investigation of many theoretical papers using the different approaches. 
The unquestionable advantage of the potential approach to the description 
of heavy mesons is its simplicity and clarity. One can calculate the masses, 
the radiative and the annihilation decay widths using a chosen potential [1-7]. 
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This allows to determine the meson wave functions in a wide range of dis- 
tances and to describe the quark interaction dynamics. The method of QCD 
sum rules has been successfully used to determine the masses of the light 
and heavy mesons [8,9]. Owing to a small relative difference of the resonance 
masses the bottomonium contribution of T-meson and ground states in other 
channels dominates only for the momenta with n > 20. The most difficulty 
of applying the sum rule method for bottomonium - is to take into account 
the relativistic corrections, i.e. the next term in the expansion in n~^. For 
this it is necessary to sum up all the terms of the type of {asn^^^)''n~^, or 
^^(asn^/^)*^ (since a^n^/^ ~ 1), which side by side with the pure relativis- 
tic effects (described by the Breit-Fermi Hamiltonian) contain the radiation 
corrections of the order of a^. This difficulty reveals, for example, in the 
calculating of masses IP-bottomonium levels. In this case the quark model 
can be used, which takes into account the relativistic effects of the quark 
interaction dynamics. 

After the discovery of the D*j{2317) [10] and Dsj{2460) [11] mesons a 
number of experimental investigations of these narrow resonances was carried 
out [10-14]. Their masses are 100 MeV smaller than those predicted by the 
quark model. A large number of theoretical papers was devoted to the study 
of structure of these mesons. The interpretation of the D*j(2317), Dsj{2A60) 
as the lowest P-wave cs-states [15, 16] is natural. However another nature of 
these states is also supposed: DX-molecule [17], tetraquark states [18, 19], 
DTT-atom [20] and others. In our study these states are considered as the cs- 
mesons. The masses of all states with the open charm and the charmonium 
states with the quantum numbers J^^ : 0~+, 1 , 0++, l+~, 1++, 2++ have 
been calculated (Table 1). The analogous computations have been carried 
out to the 6-mesons including the u, d, s, c, 6-quarks (Table 2). We have used 
the dispersion N/D method of heavy quark effective theory. 

Section II is devoted to the computation of the mass spectrum of c, h- 
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mesons with the quantum numbers J^^ = 1 , 0++, 1+ , l^'*', 2++ with 
the help of the dispersion N/D method of heavy quark effective theory (Tables 
1, 2). 

The radiative decay widths of heavy mesons are calculated in Section III 
with the help of the dispersion integration method. 

The results of the c, 6-meson characteristics computation and the prob- 
lems, which remain open, are discussed in the Conclusion. 



2. N/D method in heavy quark effective theory 



The amplitude of the N/ D method is defined as: 



AM = D(s) = 1 - B(s), (1) 



J s' - s 

(mi+m2)^ 

Here mi is the heavy quark Q{Q = c, b) mass, 1712 is the light quark q 
{q = u, d, s; qQ— meson) mass or the heavy quark Q {Q = c, b; QQ— meson; 
1712 < mi) mass, A is the cutoff parameter, p{s) is the two-particle phase 
space, which is given by formula [21-23]: 



\ (mi -|- m2) s / s 



(3) 

a, (3,5 are the coefficients, which are different for different states (Table 3), 
G{s) are smooth functions of s, which depend on the quantum numbers J^'~" 
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of the states (Table 4). The amphtude A{s) has only one singularity point 

s — (mi + m2f, therefore using the small value of ratio — ^ (m2 is the light 

mi 

quark mass) Eq.(3) can be rewritten in the form: 



-r w -r ' ' ' 

(mi + m2) 



^ . , , . I s — (mi — m^)^ . , 
We take mto account, that the expression y ~ 1 m the mid- 
dle point of the physical region (mi +m2)^ -r A. We obtain an expression for 
the function B[s), which describes accurately the heavy meson spectrum and 
has also the same structure as the explicit expression for 5-function [21-23]: 



B{s) = G{s) 



s 1 \ s - (mi + m2)^ 

a- - + i3+-d] W ^ 

(mi -|- m2j s J \ s 



A— (mi+m2)^ / s— (mi+mg)^ 



,^ Y^-y^^^^^^ ,. \ VA(A-(mi+m2)^) , 
X In V \ +i7r\ + ■ ' + (5 

' / A-(mi+m2)^ _^ l s-{mi+m2Y j [mi + m2Y ^ ' 



/AES^ (. . -^W- 



A \ ^(mi-Fm2)^ "1 ) Y A-{:mi+m2f 



Upon explicitly calculating the function B and analytically continuing it 
in s to the region (mi -|- m2)^ < s < A, we arrive at an ultimate expression 
(5) for B{s). 

By substituting the calculated function B{s) into (1), we find the roots sq 
of Eq (1). Equation (1) can be solved numerically, but it is more convenient 
to use a graphical method to find all roots sq in the complex-plane region 
under investigation. The imaginary and the real part of the expression 1 — 
i?(Res, Ims) correspond individually to a certain three-dimensional surface 
above the s plane. A section of each such surface by the s plane gives a 
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certain set of planar curves in the z—0 plane. The intersection of these sets 
of planar curves gives the set of roots Sq. Upon plotting the imaginary and the 
real part calculated by a computer individually, one finds their intersection 
and visually seeks the roots in the complex-plane region under consideration. 
The accuracy of the search may be improved by enlarging the scale of the 
graph. Considering that Sq — — iMF, we can obtain the mass M and the 
decay width F as: 

= Reso, Mr = -Imso. (6) 

The calculations show that there is only one root sq in the complex-plane 
region of interest; at a given value of A the mass and the width of the state 
corresponding to this root depend only on g that is a dimensionless constant 
from the expression for G{s). 

Four parameters were used for the calculation of the gc-meson masses: 
9qc, ^qc-i ^q-i ^c-i whcrc Qqc determines the constant Age — a parame- 

ter giving the cutoff; A^, Ac — P-wave mass shifts of the q- and c-quarks. 
The effective parameters A^ and Ac of quark mass shifts mgff = m -|- A [24] 
lead to the scale change of quark potential and allow to obtain bound states 
in P-wave. 

The values of the light quark mass (m„^d =0.385 GeV, nis =0.510 GeV) and 
the heavy quark masses {nic =1.645 GeV, nib =4.94 GeV) were taken from 
[25]. Four parameters mentioned above were determined by using the exper- 
imental data on the masses of D, D* ^ Di, Z^l^mesons. For the computation 
of cc-mesons two more parameters were used: Qc, A^ where determines 
G{s), Ac — a parameter giving the cutoff. These two parameters were de- 
termined by using the experimental data on masses of the rjc, J/^^-mesons. 
Analogous two parameters were used for the qb-mesons: Qgh, Agf,, where ggb 
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determines G{s), A^^ — a parameter giving the cutoff. Tfie parameters were 

determined by tfie experimental values of the B, B* -meson masses. 

For the c6-mesons one parameter - gcb was used, which is determined by the 

experimental value of the Sc-meson mass [26]. 

Acb is derived from Ac and A^: 

Ac6 = J(yA^+\/A^)'• 
Three parameters were used for the 66-mesons: g^, A^, A^, where deter- 
mines G{s), Aft — a parameter giving the cutoff. determines the p-wave 
mass shift of the 6-quark. These parameters are fitted on the experimental 
values of the 776(15'), T(15'), X62(l-P) -meson masses. The results of the com- 
putation are shown in Tables 1,2. It is seen from Table 1 that for D*j(2317) 
and Dsj(2460) states the calculation results on masses are 100 MeV lower 
than experimental values. It seems necessary to take into account not only 
cs, but tetraquark csnn states also. Such consideration makes possible to 
increase the masses of the D*j and Dgj states. 



3. Radiative decays of heavy mesons 

Calculations of radiative decay widths are carried out with the masses, 
which are taken from the experimental tables [26]. 

V ^-fP 

The decay width V — > 7P: {V corresponds to the notation of a vector meson, 
P - pseudoscalar meson) 



^"'^''^ ~ 24 I ^ [mv) j My |Fv(ol||Fp(0)|' 
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where a — /At: — 1/137, My is the vector meson mass, Mp is the pseu- 
doscalar meson mass. We are using the dispersion integration method [27] 



2^ f f ds' disCsdisc.,,Fv^^p{s,s' ,q^) 

"^"-^"^'^ - J ^ J V (.-M^)(.'-M|) ' 



In the center-of-mass frame [28,29]: 

Fv^,piq')r^ J dk j dyk^-[D,{s,s',,q^)Sy^,p{s,s',,q^)-eMq^)x 

where Sv^'yp{s, s', q^) is a spin factor, Gp and Gy are the vertices, which 
practically do not depend on energy, ei is the charge of quark 1, /i is its form 
factor. 



Kmax-y^^ ^ ^(.^1+1712), A - AgQ , 

where kmax is the quark momentum cutoff, 

di = — (Vl^i + k'^){ml + + A;V) , Ci ^ mj + k'^ (1 - y"^) , 
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Fv{q^)^ J dkJdyP-(^D,{s,s[,q')S'^{s,s[,q'yeMq^)x 



-1 



Fp{q^)^ J dk j dyk^-(p,{s,s\,q^)S'^{s,s[,q^)-eMQ^)x 

-1 



{s-Ml){s'^-Ml) 
The expressions for the spin factors [28, 29]: 



s' s ^ 
(s' — s)- ^ — + m| — mf) 

Sv^^p{s,s',q^) = 8<j mi 



(s'-s)2-2g2(s' + s) + 54 



s' + s — q"^ 



S$(s,s',q')^-q' 



2 2 

""2? (s' - s)2 - 2q\s' + s) + ' 



(s' + s - + 2(m? - ml)){s + s' - 2(mi - m^)) 
(s' - s)2 - 2g2(s' + s) + g4 



2 ?^ 



(12) 



-q 



(s' - s)2 - 2q'^{s' + s) + q'^ 
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/ 2 . 1 m2 — mi 

s+s —q +4mim2H x 

2 ss' 



3 \ 4g^s — (s' — s — g^)^ 
X I (s + s' - g^) ^(m2 + mi) (s + s') - ^ (m2 - mi) (s + s' - - 2(m2 + mi)^^ - 

- 2(m2 + mi)(2ss' + {ml - ml){s + s')) j] + (12) 

'^\~^'77~^ ?7^i(s+s')+"^2g^+^(m2+mi)(s+s'-g^)j 
The numerical results: 

S*(5325) ^ 7^(5279) 
Knax = 1.163 GeV, A = 39.70 GeV^ 

rB*(5325)-»7B(5279) ~ 0.8 keV 

The decay width A — > ^V: {V corresponds to the notation of a vector meson, 
A— the axial meson) 

r -^(^- (MA '\ 1 Mj + QM^M^ + M^ \FA^,vm' n 3^ 
12[ [MaJ J Ma 2M2 |Fv.(0)||F^(0)| ' ^ ^ 

where My is the vector meson mass, Ma is the axial meson mass . 

oo oo 

2^ f f ds' disCsdisCs'FA^jvis,s',q^) Q4^ 

Wto)- J - J T (s-Ml)(s'-M^) ' 

(mi+m2)2 (mi+m2)^ 

In the center-of-mass frame [28,29]: 



r^max J- 

F^^^y(g2)^ / dk f dyk^-(^D,{s,s[,q^)SA^,vis,s[,q^)-eiMq^)x 
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where Sa-^jvIs, s',q^) is the spin factor, Ga and Gy are the vertices, which 
practically do not depend on energy, ei is the charge of quark 1, /i is its form 
factor, 

FA{q')-^ J dkJdyk'-(^D,{s,s'„q')SX{s,s'„q'yeMq')x 

-1 

The expressions for the spin factors [28, 29]: 



. s' m\ + ml\fs m\ + ml\\\\ {mi + m2){s - s' + q^){s + s' - q^ + 2m\-2m\) 
^ ' 2 2 J V2 2 J j I 1 4g2s _ (s _ s' + g2)2 + 



mi-m2( s + s' -q^ , ^ ^2 i i 

^ , (17) 



s' 
where 



Zii{s,s ,q ) = 



/ 2n -(«^ + 6ss' + + ?^(2s + 2s' - g2) 



2s' 

4 

- - t-s' ^ ,r {'"■ ( (^) '-'"-^') + 

4(™lC+2(^-!!i±M)(l-^^,^ (18) 
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The numerical results: 

^,^(2460) ^7D:(2112) 

kmax = 1.139 GeV, A = 12.97 GeV^ (mi = 2030.5 MeV, ma = 565 MeV) 
^Dsj{246o)^jD*{2n2) ~ 5.73 keV. 

S ^'jV 

The decay width S — > ^V: {V corresponds to the notation of a vector meson, 
S— the scalar meson) 

2 \^ \MsJ J Ms\Fvm\Fsm' 

where My is the vector meson mass, Ms is the scalar meson mass, which is 
above threshold, i.e. Ms > (mi + m2)^. The Lehmann representation [30] is 
used here 

oo 

ds f ds' 1 disCsdisCs' Fs^-yvis , s', q^) 
~ J 'VDr{s) s'-M^ ■ 

(mi+m2)-^ (mi+mg)-^ 

The inverse propagator 

DR{m'')^ml-m^+Re{UR{ml))-UR{m^), 

where Re(llR{mR)) — n/j(m^) takes into account the finite width corrections 
of the scalar R resonance, which take into account the contribution of the 
two-particle virtual intermediate ab states to self-energy of the R resonance, 

ab 



Fs^Mq' 
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In real axis of 



/m(nij(m^)) = mrfl(m) = m ^ r{R ab, m)9{m—m.a—mb), (^23) 

ab 

where 



r{R — >• ab, m) 



^Rab 



Wnm 

is the width of the R ^ ab decay, m — mob is the invariant mass of the ab 
state, qr^^ is the couphng constant of the R resonance with the two-particle 
ab state, and 



(24) 



Pab{m^) 



1 - 



1 - 



m 



, m± = ma±mb. 



(25) 



Propagators under discussion satisfy the Lehmann representation 

1 1 7 ^^iid^) 



DR{m?) TT 



m? — w? — ie 



dm 



{rria+mi,)^ 



1 

TT 



mrji{m) 



|-Dij(m^)P(m^ — — ie) 



dm 



(26) 



(ma+mi,)2 

in the wide domain of coupling constants of the scalar R resonance with the 
two-particle ab states. 

In the one-channel case when the decay threshold is lower than the R 
resonance one has: 



Unim') = Wiim'' 



slab \ -ml) m_ mg ^ 

m+ 777.6 



167r2 I' 



771^ 



+Pa6("7^) 



i7r-|-ln 



— m_ 



— ml 



y/m^ — m'^_ 



m^ — m. 



(27) 
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at m > m+ = + m^, {nia > rrib)- 

In the center-of-mass frame: 

Fs-.^v{q^)^ dk dyk^--—x 

Jo J-i J-'r[s) 



x{D^is,s[,q')Ss-.,v{s,s[,q')-e^Mq'y^^+il^2) ), (28) 



Fs{q^) ~ / dk dye--—-x 
Jo J-1 ^r[s) 



X 



The expressions for the spin factors: 

Ss^^v{s,s',q^) = ^{m2-mi){s+ml-ml-q^){s'-{mi+m2f) - 
s 

/ 5 2 3 2 12sV+(s'-s + g2)(s + m2-m|-g2) 

X {2s'(mi(s' - s) - ms?^) + (m2 - mi){s' -s + q'^){s' - (mi + ms)^)}, (30) 



ctr/ / 2\ _ A^^i^^' ~ ^2{mi + m2){s + s' - q"^) - {mi + ■m2y{m'l - ml)) 

4g^s — (s — s' + g^)"' 

The numerical results: 

L':^(2317)^7^:(2112) 
kmax = 1.135 GeV, A = 8.94 GeV^ gn^ab ~ 1-5 GeV^. 

rD*^(2317)^7DJ{2112) ~ 2.2 keV. 

For comparison the calculation results on these radiative decays in other 
models are presented in table 5. 
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4. Conclusion 

The relativistic quark model for the study of the heavy meson spec- 
troscopy is constructed in the framework of the dispersion approach. The 
N/D method of heavy quark effective theory was used for the computation 
of the mass spectrum of the charmed and beauty mesons: the multiplets 

jpc ^ 0-+, 1 — (5-wave) and J^^ = 0++, 1++, 2++ (P-wave). The 
D*j{2317) and D<ij(2460) states attract a particular attention. The masses 
of these states are obtained to be 100 MeV smaller as compared to the ex- 
perimental data (Table 1). It seems not enough to account the contribution 
of the cs configuration only, but it is necessary to take into account a com- 
posite of tctraquark state [31] to obtain the masses close to the experimental 
data. The radiative decay widths have been calculated for these mesons using 
the experimental data on their masses. The radiative decay widths of these 
states are in good agreement with the results of other quark models (Table 
5) [32, 33]. Tabic 5 contains also the calculation results in the framework of 
the Model of QCD Sum Rules [34] and of the Vector Dominance Model [35]. 
In the present work we managed to describe all S- and P-wave multiplets of 
the charmed and beauty mesons using a small number of parameters. The 
proposed method allows to satisfy in the case of 6-mesons the qualitative 
ratio for gluon coupling constants: ggf, < gi,. The use of the dispersion N/D 
method in heavy quark effective theory allows to consider the spectroscopy 
of exotic heavy mesons and baryons taking into account the threshold singu- 
larity. 

The authors would hke to thank T. Barnes, Fl. Stancu and S. V. Chekanov 
for useful discussions. This research was supported by Russian Ministry of 
Education (Grant 2.1.1.68.26). 
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Table 1. The lowest states of the charmonium and states with open charm 
masses. 





m(0-+), GeV 


m(l~), GeV 


m(0++), GeV 


uc dc 


D 1.867 (1.867) 


D* 2.010 (2.010) 


D* 2.103 (2.352±0.05) 


sc 


Ds 1.916 (1.969) 


D* 2.105 (2.112) 


D*j 2.223 (2.317) 


cc 


ric{lS) 2.980 (2.980) 


J/ij 3.097 (3.097) 


Xco(lP) 3.393 (3.415) 




m(l+-), GcV 


m(l++), GcV 


m(2++), GcV 




uc dc 


Di 2.302 (2.420) 


L'i(2430) 2.430 (2.430) 


D* 2.460 (2.460) 




sc 


Dsj 2.350 (2.460) 


D 


,1 2.514 (2.536) 


Ds2 2.559 (2.573) 




cc 


3.726 (-) 


Xci(lP) 3.824 (3.511) 


Xc2{'i-P) 3.863 (3.556) 





The model parameters: Qqc = 1.26, A^c = 7.90, A, = 0.035 GeV, 

Ac = 0.404 GeV, = 2.92, = 5.49. 

The experimental data are given in parentheses [26]. 
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Table 2. The lowest states of the bottomonium and states with open bottom 
masses. 





m(0-+), GeV 


m(l~), GeV 


m(0++), GeV 


— — 

ub db 


B 5.279 (5.279) 


B* 5.325 (5.325) 


5.508 (-) 


sb 


5.339 (5.370) 


B* 5.419 (5.417) 


5.629 (-) 


cb 


Be 6.400 (6.4±0.39±0.13) 


6.477 (-) 


6.818 (-) 


bb 


rib{lS) 9.330 (9.300±0.040) 


T{1S) 9.460 (9.460) 


XboilP) 10.139 (9.859) 




m(l+-), GeV 


m(l++), GcV 


m(2++), GeV 




ub db 


5.534 (-) 


5.579 (-) 


5.589 (-) 


sb 


5.595 (-) 


5.664 (-) 


5.685 (-) 


cb 


6.973 (-) 


7.046 (-) 


7.073 (-) 


bb 


9.767 (-) 


Xbi(lP) 9.870 (9.893) 


Xb2{lP) 9.912 (9.912) 



The model parameters: Qqb = 3.36, Aqb = 5.70, Qcb = 2.92, gb = 4.07, 
Ab = 4.9, Afe = 0.232 GeV. 

The experimental data are given in parentheses [26]. 
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Table 3. Coefficients a, (3, 5. 



jPC 


cx 


B 




0-+ 


1 

9 


1 (mi — 

( i _1_ t ) ) ^ 2 





1— 


1 

3 


1 1 (mi — ma)^ 
6 3 (mi + 


-i(mi - ms)^ 


0++ 


i 


2 





1+- 


1 

2 


1 (mi — ma) 

2 (mi + ma)^ 





1++ 


1 
2 


1 (mi — ma)^ 

2 (mi + ma)^ 





2++ 


3 

10 


1/ 3(mi-m2)^\ 
5V 2(mi+m2)V 


-^(mi - ma)^ 
5 
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Table 4. Functions G{s). 



jPC 


G{s) 


0-+ 


8g Ag (mi + 
8 8 .s 


1— 




0++ 




1+- 


Sf? 4^' (mi + ma) 


1++ 


4^ 


2++ 


^ 

3 



Table 5. The radiative decay widths of the charmed mesons (in keV) 



Decay 


LCSR [34] 


VMD [35] 


QM [32] 


QM [33] 


[*] 


1^:^(2317) ^ 7^: 


4-6 


0.85 


1.9 


1.74 


2.2 


^,^(2460) ^ 


0.6-1.1 


1.5 


5.5 


4.66 


5.73 



The symbol [*] denotes the results of this study. 
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